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Abstract 

The leading-twist valence-quark distribution function in the pion is obtained at a 
low normalization scale of an order of the inverse average size of an instanton p c . 
The momentum dependent quark mass and the quark-pion vertex are constructed 
in the framework of the instanton liquid model, using a gauge invariant approach. 
The parameters of instanton vacuum, the effective instanton radius and quark mass, 
are related to the vacuum expectation values of the lowest dimension quark-gluon 
operators and to the pion low energy observables. An analytic expression for the 
quark distribution function in the pion for a general vertex function is derived. 
The results are QCD evolved to higher momentum-transfer values, and reasonable 
agreement with phenomenological analyzes of the data on parton distributions for 
the pion is found. 

Key words: Instanton liquid model, quark-pion coupling, pion distribution 
function and moments, renormalization group evolution 



1 Introduction 

Hadron structure functions, in terms of quark and gluon distributions specifying the frac- 
tion xp of the initial hadron momentum p carried by the active parton, play an important 
role in QCD inclusive processes. Although the evolution of parton distributions at suffi- 
ciently large virtuality Q 2 is controlled by the renormalization scale dependence of twist-2 
quark and gluon operators within QCD perturbation theory, the derivation of the parton 
distributions themselves at an initial Q 2 value from first principles still remains a chal- 
lenge. Hence, central predictions unknown in QCD are parton distributions at relatively 
low virtuality determined in a nonperturbative scheme. 
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There is some recent progress in the calculation of moments of the pion and p— meson 
parton distributions [1] within the lattice QCD (LQCD) using Wilson fermions in the 
quenched approximation, where internal quark loops are neglected. These LQCD predic- 
tions for the moments of the pion distribution function confirm the results of previous an- 
alyzes [2], being also in qualitative agreement with that extracted phenomenologically[3,4] 
from experiment [5]. However, the calculated moments are still of a relatively low accuracy. 
Besides, only a few lowest-order moments are available, while the reconstruction of the 
x-dependent distributions needs, in principle, the knowledge of all moments. Furthermore, 
the QCD sum rules calculations of parton distributions in the pion are only moderately 
successful [6], the results being justified in a limited region of the scaling variable x. 

Recently, the quark distribution function in the pion was obtained [7] in the framework of 
the Nambu - Jona - Lasinio (NJL) model [8]. These and similar studies are based on the 
assumption that the calculation of the twist-2 matrix elements, within the QCD inspired 
effective approaches, gives distributions at a low momentum scale /i ~ 1 GeV, where 
such effective theories make sense. The distributions obtained are extrapolated to higher 
experimentally accessible momentum scales using perturbative QCD, so that comparison 
with experimental data can be made. However, the problem of the NJL model is that it 
is nonrenormalizable and thus, to avoid this defect, different ad hoc assumptions about 
momentum cutoff parameters are introduced. 

The instanton model of the QCD vacuum (for recent review see , e.g., [9,10]), which gives 
the dynamical mechanism of chiral symmetry breaking and provides the solution of the 
problemfll], describes well the properties of pion [12-14] and kaon [15]. Moreover, 
it dynamically generates the momentum-dependent effective quark mass M q and quark- 
-pion vertex g nqq , and, as a consequence, provides inherently a natural ultraviolet cutoff 
parameter in the quark loop integrals through the effective instanton size p c . On these 
grounds, one may believe that the instanton vacuum framework represents an important 
advance with respect to N JL-type models. The first attempt to calculate the pion structure 
function within the instanton model has been made in [16]. More recently, important 
progress has been achieved [17,18] in calculating quark distributions in the nucleon within 
instanton-inspired approaches. 

In the present paper, based on the quark-pion dynamics in the framework of the instan- 
ton liquid model, we calculate the leading-twist valence-quark distribution in the pion 
at a low normalization point of the order of the inverse average instanton size p c . The 
calculations are performed in a gauge-invariant manner by taking into account P — exp 
factor explicitly in the definition of nonlocal quantities [19,20] and gauging the nonlocal 
quark - pion interactional, 22]. The momentum dependent quark mass and quark - pion 
vertex are constructed. The parameters of the instanton vacuum, effective instanton ra- 
dius and quark mass, are related to the vacuum expectation values (VEV) of the lowest 
dimension quark-gluon operators and to the pion low energy observables. We derive the 
quark distribution in the pion and all its moments for the general form of the effective 
quark - pion vertex function. The validity of the isospin and total momentum parton sum 
rules is ensured by the pion compositeness condition [23], and it is consistent with the 
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gauge invariant approach. As the effective instanton model is valid for values of the quark 
relative momentum up to p ~ p^ 1 ~ 0.5 — 1 GeV, the parton distributions calculated here 
are defined at this (low) normalization point p ~ p^ 1 . The results are QCD evolved to 
higher momentum-transfers, and we found reasonable agreement with phenomenological 
analyzes of the data on the pion distribution function. 

The paper is organized as follows. In Sect. 2, we briefly outline the instanton liquid model 
and introduce the quark-pion vertex. In Sects. 3 and 4, the parameters of the instanton 
vacuum model are related to the vacuum expectation values of the lowest dimension quark- 
gluon operators and to the pion low energy observables. Then, we derive the expressions 
for the moments of the pion distribution (Sect. 5) and for the ^-dependent distribution 
itself (Sect. 6), followed by the QCD evolution to higher values of the momentum transfer. 
In the last section, the results are discussed. 



2 The instanton liquid model and the quark-pion vertex 



We start with the nonlocal, chirally invariant Lagrangian of the instanton liquid model[9,10], 
which describes the soft quark fields with the soft gluon fields being integrated out. The 
corresponding action for quarks interacting through the 't Hooft vertices [11] can be ex- 
pressed in a form similar to that of the NJL model 



S inst = J d 4 x q(x)idq{x) + J d A xd 4 x'd A yd A y' K(x,x';y,y') 



1 



2N r -l 



2N r 



(qR(x')T a q L (x))(q R (y')T a q L (y)) + 



4(iV c 2 - 1) 

+ W (9*( x V<V9lO*0) (qR(y> a ^uq L {y)) +{R^ L)} . 

Here, r a = (l,rr) are the matrices for the flavor space, N c = 3 is the number of colors, 
and 

( \ 1 ± ^ 5 ( \ 
qR(L){x) = -^—q{x) 

are the quark fields with definite chirality. In the following we neglect the current quark 
mass and restrict ourselves only by the nonstrange quark sector. In eq. (1), the kernel of 
the four-quark interaction K(x, x'; y, y') characterizes the region of the non-local quark— 
(anti)quark instanton induced interaction. It is dominated by the contribution of the zero 
mode quark wave functions in the field of (anti-) instanton: 



K(x,x';y,y') — J dn(p)d 4 zdU (idip p (x — z)^ (id^ p (x' — z)^ 
(idip p (y - z)) (idipl(y' - z)) 



(2) 
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where <p p {x) is a quark zero mode profile function. In eq. (2), n(p) denotes the density of 
instantons with size p, z is the position of the instanton in the configuration space, and 
dU is the color-space phase factor. 

The spin-flavor structure of the action eq. (1) is invariant under the global axial q(x) — > 
exp (i'y^r ■ 0)q(x) and vector transformations q(x) — > exp (it • 0)q(x) and it anomalously 
violates the ?7a(1) symmetry: q(x) —> exp (ij5d)q(x). Within the instanton liquid model 
[24,13,14] it is argued that due to the long range instanton - anti-instanton interaction, 
configurations with large size instantons are strongly suppressed and the instanton density 
is sharply peaked at some finite average instanton size p c in the form n(p) = n c 5(p — 
p c ). Since the instanton liquid is assumed to be dilute, the mean separation between 
instantons is much larger than the average instanton size and the effective density n c 
is a small parameter of the approach. The values of n c and p c are estimated from the 
phenomenology of the QCD vacuum and hadron spectroscopy to be n c ~ 1 /m -4 , and 
p c changes within an interval (1.5 — 2) GeV -1 , where we put less restrictive limits on 
the range of values. It is important to note that the effective instanton size p c , which 
defines the range of nonlocality, serves as a natural cutoff parameter of the effective low 
energy model. Moreover, the coupling constants of the model, eq. (2), are also expressed 
through the fundamental parameters describing the QCD instanton vacuum, n c and p c . 
The model incorporates all attractive features of the NJL model and, at the same time, 
is free of arbitrariness in the choice of the ultraviolet cutoff procedure and physically all 
parameters are well understood. These peculiarities provide important advantages of the 
instanton model as compared to different versions of the NJL model (for a review, see, 
e.g., [8]). 

The instanton induced interaction of quarks is responsible for strong spin-dependent forces 
in hadron multiplets[25]. In particular, this force is attractive for quark-antiquark states 
with vacuum and pion quantum numbers, repulsive for the singlet part of rj', and absent 
(in the zero mode approximation) in the vector-like channels p, u, etc. If the attraction 
is sufficiently large, it can rearrange the vacuum and bind a quark and an antiquark to 
form a light (Goldstone) meson state. 

To study the formation of quark-antiquark bound states in the instanton liquid, it is 
convenient to rewrite the four-fermion term in the action eq. (1), linearizing the bilocals 
q{x)q(y) and q(x)^fq{ij) by introducing the auxiliary composite mesonQ fields M(x)[26] 
(mean field approximation) by virtue of the separability of the four-quark kernel. Then, 
we arrive at the following form of the effective nonlocal action corresponding to eq. (1): 

S = Sq + Si n t, (3) 

where Sq is the free action for quark and meson fields 



In this work we do not include explicitly the diquark part of the interaction generated by 
instantons. 
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S = J d 4 x it)q(x) + l - [a(x)(A - ml)a(x)] + ^ [n(x)(A - ml)n(x)] J , (4) 

and Sint is the quark-meson interaction part 

S int = - J d 4 Xd 4 x 1 d 4 x 2 F(x u x 2 ;nl) ' ( 5 ) 
■q(X + Xl ) E(X + Xl ;X- x 2 )[M q + g Mqq (T ■ T)M(X)} q(X - x 2 ), 

with Dirac and isospin matrices for different meson states according to (r • T) a — I • 7, 
(r ■ T) n = • t. In eq. (5), guqq is the quark-meson coupling constant and M q is the 
effective quark mass fixed in a gauge - invariant manner (see below) by the compositeness 
condition eq. (14) and the gap equation (15) in terms of the instanton density n c and 
the instanton size p c . In Eqs. (4, 5) we neglect the terms induced by tensor interaction in 
eq. (1) since they do not contribute to the scalar channels. 

To ensure the gauge invariance of the bilocal quark operators, which enters in Eqs. (3-5), 
with respect to external electromagnetic A^z) and strong A a {z) gauge fields, we include 
into eq. (5), following [21], the path-ordered Schwinger phase factors 

E(x;y) = E 1 (x;y)-E g (x;y), (6) 
S 7 (x;y) = P exp |— ieQ j dz^ A^z) J , 

E g (x;y) = P expl-tg^ j dz» A^(z)\ , 



where the charge matrix is Q = (1/3 + T/)/2, and the partial derivative is replaced 
by the covariant one = <9 M — ieA^ — igA^X 1 /2. We adopt here that the integral in the 
exponent is evaluated along a straight line with P being the path-ordering operator. The 
incorporation of a gauge - invariant interaction with gauge fields is of principal importance 
in order to treat correctly the hadron characteristics probed by external sources such as 
hadron form factors, structure functions, etc. 

The Fourier transformed gauge-invariant nonlocal vertex function F(ki, k 2 ; /j,q) describes 
the amplitude of soft transition of a pion with momentum p into a quark and an antiquark 
with momenta h\ = p + k/2 and k 2 — p — k/2, respectively. This function represents the 
full interaction vertex with all quark-gluon excitations, harder than the scale po ~ l/p c , 
strongly (exponentially) suppressed. It is defined 



F(k u k 2 ;p 2 ) = ^Q(k 1 )Q(k 2 ). (7) 
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through the quark propagator Q{p) normalized at zero to unity 

Q(p) = T^yT2 J d ^ x ex P(-iP ■ x)Q(x 2 ), (p = \p\) , (8) 



where 



Q(x 2 ) = (: q(0)E g (0, x)q(x) :)/(: g(0)g(0) :) (9) 



is the normalized instanton induced nonperturbative part of the gauge invariant quark 
propagator in configuration space. Using the explicit expressions for the instanton field 
and quark zero mode the gauge- invariant quark propagator is given by [24,19,20] 



? c 7 7 ^cos[^(arctan(^) - arctan(^))] 



o 



where R 2 = p^ + r 2 , r = \z\, t = Z4. In the derivation of these equations a reference frame 
is used, where the instanton is at the origin and x^ 1 is parallel to one of the coordinate 
axes, say /i — 4, serving as a "time" direction (i.e.,x = 0, x 4 = \x\). The propagator has 
the following expansions at small and large distances: 



Q(x 2 ) 



1 2 

1 ^r + ... as x 2 -> 0; 

2% + ... as x 2 — > 00. 

x 1 



The gauge-invariant quark propagators in configuration and momentum representation 
are plotted in Figs. 1 and 2, respectively, along with the propagators derived in neglecting 
P — exp factor in (9) and using the expressions for the quark zero mode in the singular 
and regular gauges. In the regular gauge, one has 



Qreg(% ) 2 \ 



y 2 v Vi + y 2 , 



(ii) 

y = 7T- 

2p c 



In the momentum representation, the normalized quark propagators (without P — exp 
factor) are proportional to the square of the quark zero mode in the corresponding regular 
and singular gauges: 



Qreg{p) = exp (-2p c p) (12) 
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Qsin g (p) = { Z^~ z [h{z)K X { Z ) ~ I Q (z)K Q (z)] 



PcP 

2 . 



(13) 



From Fig. 2, one can observe that in the momentum representation the shape of the 
propagator is very sensitive to the P — exp factor Q 

One of the advantages of using the gauge - invariant formalism is that the parameters 
of the model, such as the size of instantons and the effective quark mass, gain physical 
meaning. As a consequence, all other physical quantities expressed through these parame- 
ters become automatically gauge - invariant ones. Moreover, they could be compared with 
those calculated in the lattice QCD, QCD sum rules or other QCD inspired approaches. 
In contrast, when one deals with noninvariant - gauge objects there can be chosen any 
convenient gauge. It is most correct to consider the instanton vacuum field in the singular 
gauge and to construct the effective action in this specific gauge [9,10]. In the coordinate 
space in the singular gauge the instantons fall off rapidly enough to provide small overlaps 
of neighbor pseudoparticles and quasiclassical considerations are justified. But at the end 
the action has to be independent on the choice of the gauge, otherwise the form of the 
action and other observables look rather awkward. This explicitly gauge invariant form is 
suggested in Eqs. (1-5). Note also that the quark propagator, eq. (9), has a direct physi- 
cal interpretation in the heavy quark effective theory of heavy-light mesons as it describes 
the propagation of a light quark in the color field of an infinitely heavy quark [24]. 

It is important to emphasize that the meson fields entering the action eq. (3) are renor- 
malized and the field renormalization constants of composite mesons are set equal to 
zero, 



2 dU M {p 2 

9Mqq 



dp 2 



(14) 



where Um{p 2 ) is the meson field polarization operator. This condition [23] fixes the cou- 
plings of meson fields to quarks, gMqq, (see section 4) and is a consequence of the com- 
positeness of hadron states manifesting themselves as poles in the quark- (anti) quark scat- 
tering amplitude. As we will see below, it is precisely this condition supplemented by the 
gauge - invariance of the effective action, given by eq. (3), that leads to the correct parton 
isospin and momentum sum rules in the model. 



To avoid inconsistency with gauge invariance, one cannot use the treatment of the quark 
propagator in the p— representation, in factorizable form, as done in ref. [19]. 
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3 Dynamical quark mass and expectation values of quark-gluon operators 



Due to the effect of spontaneous breaking of the chiral symmetry, the momentum de- 
pendent quark mass M q (k) is dynamically developed. It obeys the well-known gap equa- 
tion [9, 14] Q 

d'k M*(k) n c 

(15) 



(2tt) 4 k 2 + M 2 (k) AN ( 
the solution of which has the form 



M q (k) = M q Q{k). (16) 



Given the dynamical mass, the values of the quark condensate, 

A , T f d 4 k MJk) 



and the average quark virtuality in the vacuum [27], 

>2 _ (: gD 2 q :) _ AN C f d*k fc2 M q (k) 



q 



(:qq:) (qq) J (2vr) 4 k 2 + M 2 (k\ 



can be found. The average quark virtuality defines the derivative of the quark conden- 
sate and thus nonlocal property of it. One of the main suggestion of the QCD sum rule 
method [28] was that the local quark and gluon condensates dominate in the light hadron 
physics and introduction of higher dimensional corrections or even nonlocal condensates 
themselves [27] have not to change the standard results too much. Thus, at least for con- 
sistency of local and nonlocal QCD sum rules, the derivative (virtuality) value has to be 
relatively small. Phenomenologically, there is rather fine QCD sum rule analysis of this 
value [29] where it has to be defined as X q « 0.4 ± 0.2 GeV 2 . The LQCD calculation 
yields X q = 0.55 ± 0.05 GeV 2 [30]. Certainly, there is corrections from direct instantons 
to the QCD sum rule result, but they have not to change the result drastically. It would 
also be urgent if the LQCD estimation could be confirmed by new calculations. The ratio 
7] = M q (\ q )/\ 2 characterizes the diluteness of the instanton liquid vacuum. The smallness 
of rj means that the dynamically generated quark mass is not big enough essentially to 
modify the instanton vacuum. 



3 Here and in the following, all Feynman diagrams are calculated in the Euclidean space 
(k 2 = —k"fi) where the instanton induced form factor is defined and rapidly decreases, so that 
no ultraviolet divergences arise. At the very end we simply rotate back to the Minkowski space. 
One can verify that the numerical dependence of the results on the pion mass and the current 
quark mass is negligible and can be dispensed with in the following considerations. 
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For the moment, it is instructive to consider equations Eqs. (15 - 18), neglecting the 
term M 2 (k) compared to k 2 in the denominator of the integrands. This approximation is 
justified in the dilute liquid regime where (k 2 ) = X 2 » M 2 (X q ). The observed accuracy 
of such procedure is better than 20 — 30% if rj < 1. Then, from Eqs. (17) and (18), by 
using the explicit forms given in Eqs. (10), (8), (16), we have 

N c M q 2 2 

{qq) = " q = p7- (19) 



The first relation, put in the form M q = —(2ir 2 /N c )p c 2 (qq), coincides with the result 
obtained in ref. [31], where the effective quark mass has been defined in a system of 
small size instantons interacting with long wave vacuum fields. The coefficient in this 
relation is equal to the normalization factor of the momentum representation of the quark 
propagator. It turns out that this factor, which is equal to (2n) 2 , as seen in Eq.(8), is 
the same which appears in the gauge-invariant propagator and also in the singular gauge 
propagator (without Pexp). 

The second relation in eq. (19) has recently been obtained in ref. [20], where nonlocal 
properties of the quark condensate are studied within the instanton model. The same 
result was also obtained in ref. [32] from direct calculations of the local mixed quark- 
gluon condensate in the framework of the Diakonov - Petrov model: 



\ 2 



X a 

(: q{ig(T^Gl v —)q :) 

Vm 7 ) 



(20) 



It is clear, from the expressions for the average quark virtuality, that the range of the 
quark - antiquark interaction is characterized by the effective size p c of the instanton 
fluctuations in the QCD vacuum. The natural gauge - invariant definition for the average 
quark virtuality, eq. (18) (and also that in eq. (17) for the quark condensate), with M q (k) 
defined in Eqs. (8) and (16), is valid only if the zero mode solution in eq. (8) is written 
in the gauge-invariant way. If we substituted its expression in the singular gauge (in 
neglecting Pexp factors) in eq. (18), we would obtain A^ = 9/(2p^), with a coefficient far 
from the correct one. 

Inverting the relations eq. (19), we express the parameters of the instanton vacuum 
model in terms of the fundamental parameters of QCD vacuum p 2 = 2/X 2 , M q = 
— (ATr 2 /N c )((qq)/\ 2 ). If we expected "standard" values for the quark condensate, ( qq) 
-(230 MeV) 3 (see, e.g., [10]), and for the average quark virtuality, X 2 q ps 0.5 GeV 2 [29,30], 
we would obtain p^ 1 ps 0.5 GeV and M q ps 0.32 GeV. As we will see in the next section, 
the joined analysis of the vacuum and pion properties confirms this estimation. 

Within the dilute liquid approximation, the gap equation, (15), leads to 
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n c = 



N C M 2 \ 2 



4tt 2 



*-I n , with I n = J duuQ 2 {u/p c ) « 0.61, (21) 



where the constant I n is independent of p c . There are other different useful combinations 
relating vacuum parameters with each other. For example, 

<«> = ~hk = -^7T^ (22) 



These relations have the same parametric dependence as in the estimations obtained in 
[24,10] but with different coefficients. The first one expresses the quark condensate in 
terms of the effective single instanton contribution times the density of instantons. The 
reason for the difference in the coefficients is that in [24], where it looks as (qq) = —n c /M q , 
the expressions were obtained from the instanton formula in the gas approximation by ad 
hoc replacing the current quark mass by the effective quark mass. In contrast, in deriving 
eq. (22) this replacing procedure is fixed by the gap equation eq. (15), with a definite 
coefficient mainly defined by the slope of the form factor M q (k). The second relation 
represents a self-consistent value of the quark condensate in the instanton vacuum model 
(cf. [10]). Further, since the instanton contribution to the value of the gluon condensate 

is given by: (— G 2 )\i nst = 8n c , it can be expressed through the quark condensate and the 

7T 

average quark virtuality 

<-G 2 >U = ^iM! < 0.019 GeV\ (23) 



The "standard" value of the gluon condensate estimated in the original work, in ref. [28], 
was ((a s /ir)G 2 ) ~ 0.012 GeV 4 . The latest reanalysis [33] of the QCD sum rules for heavy 
and light mesons and also recent lattice results [34] provide values which are twice or even 
larger than the "standard" one. 



4 Pion low energy observables 



Let us now consider the low-energy observables of the pion. The pion - quark coupling 
constant is determined by the compositeness condition, eq. (14), with the pion mass 
operator being 

U 7T (p 2 )=N c J -0- 4 F 2 (k,k + p-,p 2 o )Tr{ l5 S(k + ph,S(k)}, (24) 

where the normalized nonlocal vertex is given in (7) and the quark Green function is 
S(k) = [M q Q(k) + k}- 1 with Q(k) defined in (8). 
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From the definition eq. (14) we derive the expression for the pion - quark coupling constant 

9nqq 



9nqq ~ t\t T I 2 \ ' (^) 



In the case of a massless pion, the integral I gn reduces to 



V(o) 



dk k 3 Q(k) 
D 2 {k) 



kQ'(k) (kQ'(k)Y 
4Q(k) \4Q(k)J 



(26) 



where 



Q'(k) = -^Q(k), D{k) = MlQ\k) + k 2 . 



(27) 



The expression for g^qq given in Eqs. (25) and (26) agrees with that derived in [14]. 

To fix the parameters in the instanton model, we consider the low energy decay constants 
of the pion. As it has recently been shown in [22], in the presence of nonlocal separable 
interaction the axial current conserved in the chiral limit can be constructed from the 
action eq. (1) by using a Noether - like method.^ The full current is the sum of local, 

J5(L)(^) = ^(x)7 M 75r a g(x), (28) 



and nonlocal, 



jh(ni)( x ) = / 11^* K(x 1 ,x 2 ,x 3 ,X4){A 1 (x 1 ,X2,x 3 ,x 4 )q(x 1 )i'y 5 T a q(x 3 ,)q(x2)q(x i ) + 
J i=i 

+A 2 {x 1 ,x 2 , x 3 , x A )ie ahc q(xi)T c q(x 3 )q(x 2 )i^^ h q{x i )} (29) 

pieces. The coefficients A\ and A 2 are derived in ref. [22] and in principle depend on 
the gauge fixing procedure. Fortunately, there is no path dependence for the longitudinal 
components of the current, and thus, the decay constants considered below are well- 
defined. 

The axial and vector currents in different isospin states have a similar structure [22]. As 
a result, various Ward identities which follow from (partial) current conservation and the 

4 One of us (A.E.D.) thanks M.C. Birse for discussion of the problem of current conservation 
in the nonlocal models. 
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low-energy theorems are satisfied. In particular, the Goldberger - Treiman relation for the 
quark-pion coupling constant has the usual form 



dirqq = "p ( 30 ) 
Jit 



and the 7r° — > 77 decay constant 



^ 77 "47T 2 / 7r (31) 



is consistent with the requirement of axial anomaly. 

We fix the model parameters to give the pion weak decay constant f n , within 1% of 
accuracy. In Table 1, we present the results. For the two model parameters, M q and p c , 
we show the predictions for the quark-pion coupling g nqq , the quark condensate (qq), the 
average vacuum quark virtuality \ 2 and the instanton density n c . 

Table 1. The values of the low energy vacuum and pion observables discussed in the 

text. 



(GeV) 


Pc 

(GeV- 1 ) 


U 
(MeV) 


9nqq 


1(^)1 1/3 
(MeV) 


(GeV 2 ) 


n c 
(fm" 4 ) 


0.18 


1.0 


92 


3.7 


294 


2.2 


1.56 


0.22 


1.5 


92 


5.5 


233 


1.0 


0.86 


0.23 


1.7 


93 


6.2 


215 


0.83 


0.70 


0.26 


2.0 


91 


7.7 


197 


0.63 


0.57 



As shown in Table 1, the values of the parameters M q and p c that reproduce the lowest 
dimension VEV with an accuracy of an order of 30 % are in the "window" M q = (220—260) 
MeV, p c = (1.5 — 2.0) GeV -1 . In the following we use the typical set of parameters: 

M q = 230MeV, p c = 1.7 GeV- 1 . (32) 



The diluteness condition rj = M 2 (A g )/A 2 << 1 is well satisfied within the whole "window". 

The momentum dependence of the vertices in the numerators of the integrands (which are 
defining physical quantities) is important because it provides the ultraviolet regularization. 
Also, due to momentum dependence of the vertices, the measure in the integrals looks like 
product of some powers of k 2 and the function Q(k). This measure has maximum at k 2 of 
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order 1/ p 2 c and, at small momenta, the momentum dependent quark mass in denominators 
can be substituted by an effective constant mass parameter m q pa M q (k ~ p^T 1 )- With the 
form of momentum distribution shown in Fig. 2, it approximately equals to the mass at 
zero, M q (0). This mass parameter m q has to be identified with the standard constituent 
quark mass. Corresponding to this substitution, we redefine the function D(k) given in 
eq. (27) as D(k) = m 2 + k 2 . The choice of the mass parameter, 

m q = 230 MeV, (33) 

well reproduces the integrals defining the VEV given by Eqs. (15 - 18). This constant - 
mass approximation is often used in practice with the quark mass in the region 250 — 350 
MeV (see, f.i. [18,35,36]). 

The model parameters and predictions for vacuum and pion observables are obtained 
within a set of approximations. We are working in the chiral limit of zero current quark 
mass. Further, within the zero mode approximation small contributions coming from 
vector mesons are neglected. Also only the lowest two - quark Fock intermediate state 
in the pion is taken into account, which corresponds to the quenched approximation. We 
regard that all these factors can change a little the numbers in Table 1, but the qualitative 
results discussed are not greatly influent. 



5 Moments of the quark distribution function 

The standard QCD analysis based on the Operator Product Expansion (OPE) relates 
moments of parton distributions at a given scale to the hadronic matrix elements of 
local twist-2 operators. This formalism is employed to separate the hard and soft parts 
of the forward scattering amplitude. Within the OPE, the hard part is calculable within 
perturbation theory in the form of Wilson coefficients. The soft part is represented by a set 
of local operators classified by the twist. Their matrix elements accumulate information 
on the nonperturbative structure of the QCD vacuum. 

The twist-2 gauge - invariant non-singlet local quarkF] operators with flavor j are defined 
by 

O j w . m = t N q,{ llll D, 2 ..D llN } s q v (34) 

where = — igA^Ta is the covariant derivative and the symbol {...}s means the 
traceless and symmetric part of the tensor. The matrix elements A 3 N of the local operators 

5 As in ref. [18], we will neglect gluon operators, justified by the smallness of the diluteness 
parameter r\. 
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J N between pion states |7r(p)) with momentum p, renormalized at the normalization scale 
fi, are defined by 

= Y (n(p)\q ] n(n I/ Dn N q Mp))\», (35) 



where n v is a light-like vector, with n 2 = and (np) = 1, introduced to select the symmet- 
ric traceless part of the operator J N , eq. (34). Let us now define the quark distribution 
for the j—th flavor in terms of its moments, viz.: 

l 

A j N {fi 2 ) = Jdx x N - 1 q j (x,n 2 ). (36) 
o 



The x— variable is the fraction of the longitudinal pion momentum carried by a quark in 
the infinite-momentum frame. The fi 2 dependence of A J N is known exactly from the solu- 
tion of the perturbative QCD evolution equations, while the nonperturbative dynamical 
model provides the initial input for this evolution. These initial values of the moments are 
calculated here in the instanton model, which specifies a low momentum transfer value 
related to the scale /Xq oc 1/ 'pi. 

The N— th moment A v ^\jj^) in the instanton model, from Eqs. (34 - 36), can be expressed 
as (see Fig. 3) 



d A k 

W) 4 



(37) 



{F 2 (k, k + p; f4)Tr[ l5 S(k + p)^ (k + p), 2 ..(k + p) flN } s S(k + p) l5 S(k)]- 
_ 2 fF*(k, k+p; tih Tr[lbS{k +p){(k +pMk +p) ^ (k +p) IXN } sl5 S(k)\ 



d(k + p) 2 

Here, due to the nonlocal character of the interaction Eqs. (5-6), the additional term with 
a derivative ensures gauge invariance of the approach and enables us to satisfy the isospin 
and momentum conservation sum rules. Indeed, taking into account the compositeness 
condition eq. (14) we get for the first two moments 



dU n (p 2 



dp 2 



= 1, <(^) = ^ 



dU n (p 2 



dp 2 



1 

2' 



(38) 



The results for the first two moments, eq. (38), manifest the (normalization) isospin and 
momentum sum rules for the valence-quark distribution function as 

l l 
J dx q val {x,fi 2 ) = 1 and J dx x[q val (x , fi 2 ) + q" al (x,n 2 )} = 1, (39) 
o o 
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where q val (x) is the antiquark distribution. The fact that, at the low momentum scale 
/j,q, the whole momentum in the pion is carried off by the valence quarks is due to the 
quenched approximation used, when only valence quark- antiquark intermediate states are 
included and all intrinsic quark-gluon sea states are neglected. 

Thus, in the quenched approximation the dynamical information contained in the first 
two moments is strongly restricted by the symmetries and kinematics, and as a result, it 
is model independent. The nontrivial dynamics is contained in the moments with N > 2. 
The general structure of the moments of the structure function (SF), from eq. (37), can 
be written in the form 

-«^f'^(V)™' *=1,2,... (40) 
with the coefficients Jf F given by 

where the vertex terms with derivatives, like that appearing in eq. (26) for I gn (0), are 
denoted by dots. In eq. (40), the square brackets [...] mean the integer part of the number, 

and are the binomial coefficients. 

It is instructive to consider two extreme cases, depending on the physics under consid- 
eration. If the QCD vacuum were a very dense medium, r] >> 1, then jf F = for all 
i except i = 0. As a result, it leads to the set of moments A N = 1/2^ N ~^ for all n and 
to a quark distribution which has the form of a delta function: q(x) = S(x — 1/2). This 
extreme case corresponds to the heavy quark limit, and the coefficients Jf F represent 
consequent corrections in inverse powers of the heavy quark mass: ~ ((k 2 )/m 2 y. In the 
opposite extreme case of a very dilute vacuum rj << 1 one gets jf F = 1 for all % and 
Ajy = 1/-/V for the moments. This extreme case corresponds to the momentum indepen- 
dent quark mass and provides flat quark distribution q(x) = 1. Moreover, the first term in 
eq. (41) dominates over the terms indicated by dots, since the latter are small of an order 
of 0(p c m q ). A realistic situation seems to be somewhere in-between these two extremes. 
Note that the role of pion mass is negligible, but the interplay of the effective quark mass 
and the slope of the nonlocality in Q(k) has an important effect. 



6 Quark distribution function and its QCD evolution 

Let us now turn our attention to the quark distribution itself. This distribution for the 
pion with 4-momentum p is given by (see a graphical representation in Fig. 3) 
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q(x; nl)p» = 2N c g 



d 4 k 



lm I J^yA* - 1 - ■ n)] 



■Tr | 75 5(A; + p) F*(k, k + p; ^S(k + p) - 2 (^|^^) (* + P), 



(42) 

75^)1, 



where q(x) = u(x) va i = d(x) va i for ir~ . Here, we arrive at the quark distribution defined 
in a similar manner as that used in [37,38]. The S[x — 1 — (k ■ n)] — function appearing 
in eq. (42) represents the effective vertex related to the composite operator 3 N given by 
eq. (34). It accumulates information about all moments of the distribution function (37) 
and is related to them by a Mellin transformation if the light-cone vector n is normalized 
by (pn) = 1. The light-cone vector n serves to project out in Eqs. (37) and (42) symmetric 
traceless tensors. It can be easily shown that the first moments of q(x) will reproduce the 
parton sum rules eq. (38). 

To calculate the A;— integral in eq. (42), we use a— representation for the propagators^, 



k 2 + w? 



da exp — a(k 2 + m 2 ) 



(43) 



and for the vertex 8— function, 



1 °° 

5[x — (k ■ n)] — — / da exp [ia(x — k ■ n)]. 

2,71 J 



(44) 



Then, a direct calculation from eq. (42) provides the result for the quark distribution, 



which in the massless case (mi 



0) is reduced to 



N c 9l qq 



2tt 2 





' m 



J j dv x dv 2 F{v^)F{vi) exp (- 

Q(xv 2 > xui) + (x <-> x) \ . 



(45) 



Ei I — - + x exp 

XVi) 



XVi 



In the above equation, x — 1 — x, Ei(z) is the integral exponential, and F(u) is the 
correlation function related to the vertex function Q(k) by the Laplace transformation. 
The vertex Q(p), in the essential region of p (0 < p < 4:/p c ) is approximated by 



Q(p) = 4.5 exp (-1.9p c p) - 3.5 exp (-3.6p c p), 



(46) 



For details, see ref. [38]. 
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which leads to the Laplace transform, 



Fly) 



8.55exp(-0.9p^z/) - 12.6 exp(-3.24p^z/) 



(47) 



Let us stress that the expressions eq. (40) for the moments and eq. (45) for the valence- 
quark distribution in the pion are universal ones and valid for any shape of the functions 
Q(k) and F(u), which in turn are determined by the concrete model of the quark-pion 
dynamics. 

The quark distribution q(x, /Jq) and the momentum distribution (structure function) 
xq(x,[iQ) are shown graphically in Fig. 4. We have to note that the shape of the dis- 
tribution is quite stable with respect to changes of the instanton model parameters, if 
they are fixed to reproduce the pion low energy properties. We remind that we are com- 
puting only the leading-twist distributions at a low normalization point /io ~ p^ 1 rather 
than the full structure functions which contain also higher-twist corrections. The latter 
may be large at low q 2 . We have also to note that these results differ strongly from those 
obtained in calculations with the NJL model [7] which yield distributions that are rather 
consistent with the strict chiral limit q(x, /ify ~ 1. 

The computed distributions are then used as initial conditions for the perturbative evolu- 
tion to higher values of Q 2 , where the power corrections are expected to be suppressed, so 
that one can compare them with the available experimental data. Actually, we compare 
our theoretical predictions with the phenomenological analysis by Sutton et al [3] of the 
data taken from Drell - Yan and prompt photon experiments performed by the groups 
NA10 (CERN) and E615 (Fermilab) [5]. 

The form of the evolved distribution q(x, Ql) at the momentum scale Ql = 4 GeV 2 is 
reconstructed from its moments evolved to this scale in the leading order (LO) and next-to 
-leading order (NLO) perturbative QCD in the MS scheme by using the first six Jacobi 
polynomials. To this goal we use the well-known expressions [39] for the perturbatively 

calculable coefficient function of the process = + — - C^ and the anomalous 

dimensions 7( n ) calculated up to LO and NLO. Thus, the final result for the moments 
obtained from the factorization procedure is 

i 

An(Q 2 ) = £Cf (QV)OfV) = fdx x N q(x,Q 2 ). (48) 
< o 



In performing the evolution analysis we choose a low momentum scale /i = 0.19 ± 
0.05 GeV 2 , and a set for the QCD scale parameter Aq^ d = 0.19 GeV in order to be 
consistent with [3]. The resulting distribution q(x, Qq) is shown in Fig. 5 together with 
the phenomenological curve derived from the data in [3]. The initial distribution function 
at the low-momentum scale /Xq is also shown for comparison. 
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The values of the first moments of the pion quark distribution at Qq = 4 GeV 2 calculated 
in LO and NLO are shown in Table 2. The error bars quoted in Table 2 for our calculations 
are due to accepted uncertainty in the choice of the initial scale of evolution /x . These 
values should be compared with those obtained from the phenomenological analysis [3] and 
from LQCD simulations [1]. In Table 2 we also include the moments of quark distribution 
in the pion obtained from the parametrization [4]. 



Table 2. The values of the first moments at Ql = 4 GeV 2 . 





LO 

(this 


NLO 
calculations) 


LQCD[1] 


Exp. fit [3] 


Exp. fit [4] 


MQl) 


0.321 ±0.02 


0.293 ±0.03 


0.279 ±0.083 


0.230 ±0.01 


0.193 ±0.01 


MQl) 


0.150 ±0.02 


0.112 ±0.025 


0.107 ±0.035 


0.101 ±0.005 


0.083 ±0.005 


MQl) 


0.083 ±0.015 


0.057 ±0.025 


0.048 ± 0.020 


0.057 ±0.005 


0.046 ± 0.005 



Let us finally discuss the uncertainties of the QCD evolution from the low momentum 
scale hq. As we see from Table 2, the difference of the LO and NLO results is in the range 
of 30%. It turns out that the use of a larger initial evolution scale, say /Iq > 0.3 GeV 2 , 
gives a rather good convergence with deviations less than 10%, whereas in the opposite 
case, i.e., for scales smaller than about 0.1 GeV 2 the deviations increase and perturbative 
evolution loses any sense. This behavior has also been observed in analyzes within the 
NJL model [7]. 

The comparison shows that our calculations, in particular in NLO, are consistent with the 
phenomenological analysis of [3] and fairly close to the LQCD results. Both theoretical 
approaches (LQCD and the instanton model) predict moment values systematically larger 
than the phenomenological one. One of the reasons for this disagreement may be traced to 
the quenched approximation which does not take into account any sea-quark contributions 
at the initial scale, attributing in this way the whole pion momentum to the valence quarks. 
Indeed, the origin of the A 2 moment at the initial scale (in the quenched approximation) 
and its subsequent evolution is purely kinematic and does not depend on the details of 
the model. In principle, one could match the valence momentum fraction derived in our 
calculation with that determined in [3] by shifting the initial value /i 2 , down to 0.01 GeV 2 
(see, for instance, [7] a). However, to start a perturbative evolution from this very low 
scale is formally incorrect and technically amounts to a rather unstable procedure. In our 
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opinion, it is more realistic to expect that including in our analysis contributions from the 
sea, which may have momentum fractions of an order of 10%, the agreement between the 
theoretical predictions and the phenomenological analysis can be considerably improved. 
In addition, the effect of nonperturbative evolution [40] from an initial scale //q up to 
Q 2 ~ 1 GeV 2 could be important. 



7 Results and discussions 

In summary, we have presented theoretical predictions for the valence-quark distribution 
function, eq. (45), and its moments, eq. (40), for the pion. The calculations are based 
on the instanton model of the QCD vacuum as a candidate treatment of nonperturbative 
dynamics, expressing the observable hadron properties in terms of fundamental character- 
istics of the vacuum state. We found that the instanton model describes well the vacuum 
expectation values of the lowest- dimension quark-gluon operators and the pion low-energy 
observables. To obtain these results, we have used gauge invariant forms for the dynami- 
cally generated quark masses and quark pion vertex. Thus, we are led to express the form 
of the pion quark distribution function in terms of the effective instanton size p c , and the 
quark-mass parameter m q . The pion quark distribution function extracted corresponds to 
a low normalization scale, where the effective instanton approach is justified. It is shown 
that the validity of parton sum rules for the isospin and total momentum distribution 
is a consequence of the compositeness condition and the strict implementation of gauge 
invariance. We have used techniques to derive these results which constitute a complemen- 
tary approach to lattice simulations and to phenomenological fits to experimental data. 
Using this distribution function as an input, we obtained the quark distribution function 
in the pion via standard perturbative evolution to higher momentum values, accessible 
by experiment. A reasonable agreement with the data was found. In fact, the calculations 
are performed in the quenched approximation, where the effect of intrinsic quark-gluon 
sea is neglected. We expect that the effects of the intrinsic quark component of the pion 
wave function and the nonperturbative evolution at intermediate energy scale provide a 
better agreement between theoretical predictions and phenomenological analysis. 
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x /p 



Fig. 1. Configuration space representation of the normalized instanton induced nonperturba- 
tive part of the gauge-invariant quark propagator, eq.(10) (solid line); and the corresponding 
propagators derived without P — exp factor in the singular (short-dashed) and regular, eq.(ll), 
(long-dashed) gauges. 
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Fig. 2. Normalized momentum space representation for the same propagators given in Fig.l, 
corresponding to eqs. (8) (solid line), (12) (long-dashed) and eq.(13) (short-dashed). 
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k — p 



Fig. 3. Graphical representation of the operator product expansion. The left hand side of 
this diagram is the imaginary part (.Discontinuity) of the forward scattering amplitude. Within 
OPE it is represented by the convolution of the Wilson coefficient function Cns(Q 2 ) of a "hard" 
parton subprocess (upper block of the right diagram) and the "soft" parton distribution function 
q(x,Q 2 ) (lower block of the right diagram). The constituent quark and pion are depicted by 
solid and double lines, respectively. The wavy line denotes the virtual photons. Ons is the local 
operator and cross on the quark line correspond to 5{x — (kn)). 
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Fig. 4. The valence-quark distribution function (DF) q(x; /Xq) (dashed line) and the quark 
momentum distribution function (MDF) xq(x; (solid line) for the pion as a function of the 
longitudinal momentum fraction x at the low momentum scale /j,q = 0.19 GeV 2 and density 
parameter p c m q = 0.39. 
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Fig. 5. The quark momentum distribution function xq(x; Qq) (long dashed line) for the pion as 
a function of the variable x evolved to the momentum scale Qq = 4 GeV 2 (LO approximation) , 
using p c m q = 0.39 for the density parameter. The solid line denotes the phenomenological curve 
[3] on the same scale Qq, extracted from the data. The short - dashed line shows the same 
distribution on the low momentum scale /Xq = 0.19 GeV 2 . 
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